The shape of disorder broadened Landau subbands in graphene 
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Density of states (DOS) of graphene under a high uniform magnetic field and white-noise random 
potential is numerically calculated. The disorder broadened zero-energy Landau band has a Gaussian 
shape whose width is proportional to the random potential variance and the square root of magnetic 
field. Wegner-type calculation is used to justify the results. 
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Graphene is a two-dimensional (2D) system that has 
attracted a lot of attentions in recent years [Ull] because 
of the distinct and similar electron properties between 
graphene and conventional 2D electron gases such as the 
inversion layers of semiconductor heterostructures. Dif- 
ferent from the usual 2D electrons in solids that obey 
Schrodinger equation, the governing equation of electrons 
in undoped graphene is the relativistic massless Dirac 
equation ,3| due to the linear dispersion near the vicin- 
ity of the two nonequivalent corners K and K' of the first 
Brillouin zone. Similar to the 2D Schrodinger electrons, 
2D Dirac electron states in a perpendicular uniform mag- 
netic field B also form highly degenerated discrete Lan- 
dau levels (LLs). Instead of equal spaced LLs for the 
Schrodinger electrons, the LLs of 2D Dirac electrons are 
unevenly distributed 0, Q according to 



E n = ±j2ehv 2 F nB 



(n = 0,1,2...), 



(1) 



where vf is the Fermi velocity, e is electron charge and 
h is the Planck constant. The existence of zero-energy 
Level for n = is a direct consequence of charge-hole 
symmetry. In the presence of disorders, LLs broaden 
into Landau subbands, gives raise to the quantum Hall 
effects Q. The distribution of electron levels inside the 
Landau bands is a subject of fundamental interest be- 
cause it is important to physics quantities sensitive to the 
density of states (DOS). For 2D Schrodinger electrons, 
the subject have been studied both theoretically and ex- 
perimentally for a long time@,S,[l[iS[ll[il[iE[il[ii. 
Based on self-consistent Born approximation (SCBA), 
the shape of DOS was argued to be the well-known semi- 
circle formQ. Later, Wegner predicted that the shape 
of the lowest LL in the strong-field limit and with white- 
noise random potential^ is Gaussian-like. The super- 
symmetry method reaches the similar conclusions (1 Ot 1 1 2j | . 
However, the shape of disorder broadened Landau bands 
for the Dirac electrons is less well studied although it 
is important for the issues of diagonal and non-diagonal 
conductivity [3]. In this paper, we show that the shape of 
zero-energy Landau band is best described by Gaussian 



function. The broadening width T of the band depends 
on both the magnetic field and disorder randomness. In 
the high magnetic field, it is proportional to the square- 
root of the field. At a fixed field, the broadening width 
is approximately linear in the square-root of disorder po- 
tential variance. Interestingly enough, these results are 
similar to the predictions of Wegner [8| on Schrodinger 
electrons. To demonstrate the importance of the shape 
of Landau subband in electron transport in high mag- 
netic field, the magneto-conductivity is calculated. 

The it— electrons of graphene in a perpendicular mag- 
netic field can be modeled by a tight-binding Hamiltonian 
on a honeycomb lattice of two sites per unit cell, 



\i>< j\ + h.c.)+J2 £ * 



i >< i\ 



(2) 



where t — —2.7 'eV is the hopping energy between the 
nearest neighboring atoms. To mimic randomness on 
Dirac electrons, the on-site energy is assumed to be 
random with Gaussian white-noise distribution of zero 
mean < Si >= and finite variance < SiEj >— Sijg 2 . 
g measures randomness. The magnetic field B is in- 
troduced by means of Pcricrls' substitution in hopping 
parameter 13, [HI as t — > te* 3 ', where 4>i — 2ne/h J J A- 
dl. 

In order to obtain numerically an accurate averaged 
DOS of the above Hamiltonian, the system size should 
be large enough so that the Lanczos recursive method (l9| 
is used. The averaged DOS can be evaluated by p{E) = 
< ip\Im E _lj +iv \ip >, where v is an infinitesimal positive 
number. The bar denotes the ensemble average. When 
disorder is present, the system does not have transla- 
tional symmetry and all sites are no longer equivalent. 
Estimating the average by using the " site-basis" is there- 
fore inefficient. Inspired by the work of TriozonJlO], we 
use the random phase state [2lj. Thus we only need to 
consider several random phase states to get a good es- 
timate of the average DOS. In this approach, a small 
artificial cut-off energy (e = lmeV) is introduced to sim- 
ulate the infinitesimal imaginary energy v[22, 23]. This 
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FIG. 1: (Color online) (a) DOS in a magnetic field of 10T 
(solid line) and DOS without magnetic field (dashed line) for 
various g. (b) DOS at B = 10T for 3 = 1,2, 3, 4eV. 



artificial parameter will lead to a small width of LLs in 
clean graphcnc. 

Both on-site energy randomness and the magnetic field 
can affect DOS. There are two important length scales. 
Disorder causes electron scattering, leading to an elec- 
tron mean- free path I. The magnetic field introduces a 
magnetic length l c — y/h/eB. Electronic states in a Lan- 
dau level (LL) can be viewed as cyclotron motions of an 
electron around orbits of radius l c centered at different 
location. Scattering can destroy and modify these cy- 
clotron motions, leading to Landau subband broadening. 
Thus, one can expect different DOS at two limits, l c 3> I 
and l c <C I- Fig. la is the DOS for l c ^s> I. Indeed, our 
numerical results show that the DOSs at B = 10T and 
strong disorder (in the sense of l c 3> I) are almost the 
same as DOS at zero field and finite g, showing no obvi- 
ous magnetic field dependence of the Landau subbands 
in the limit. Fig. lb shows DOS of different g at a fixed 
field B = 10T at which l c 3> Z is still satisfied. The van 
Hove singularities is less singular as g increases whereas 
the DOS near the Dirac point (E = 0) increases with g. 
These results agree with prior analytical results 24]. In 
the opposite limit (high magnetic field), the strength of 
random potential becomes relative weak such that Dirac 
electrons can complete their cyclotron motion before en- 
countering scatters. For g = O.leV, we observed the os- 
cillation of the DOS near the Dirac point as shown in 
Fig. 1(a), which reflects the appearance of LLs in high 
magnetic field. In what follows, we focus on the Landau 
subband broadening in high magnetic field. 

In the limit of l c <C I, DOS shows a very different be- 
havior. Fig. 2a is the DOS of graphene without disorder 
(g = 0) and in a uniform magnetic field (B = 10. 0T). 
<5 — like spectrum distributed according to Eq. (1) is ob- 
served. The seemingly non-zero width is due to the fi- 
nite size effect and cut-off energy lmeV introduced in 
our numerical calculation. Furthermore, these finite-size 
broadened spectrum can be fitted well by Lorentizian dis- 
tribution of 2meV width (from the cut-off energy). This 
artificial broadening may also be viewed to be a true 
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FIG. 2: (Color online) DOS at B = 10T. (a) Clean graphene 
g = 0; (b) disordered graphene with g = 0.05 0.2, 0.3 0.4eV; 
(c-e) circles are numerical results of DOS of the zero-energy 
Landau band for g = 0.05eV (c), #0.2eV (d), and 0.4eV (e). 
The solid (red) lines are the Gaussian fitting of the form 



= exp{- 



Ty'ir/2 

pure graphene. 



-). The dotted lines are the zero-energy LL of 



broadening accounting the electron-phonon or electron- 
electron interactions in reality. 

Fig. 2b is the DOS of Landau subbands at B = 10T for 
various g. The LLs degeneracy are lifted in the presence 
of disorder, LLs are broadened into Landau subbands. 
The broadening width (r) of zero-energy subband in- 
creases from 3.0meV to ll.OmeV when g changes from 
0.05eV to 0.6eV. When the strength of the white-noise 
potential is very weak (g < O.leV), the broadening width 
of each Landau band increases a little. As shown in Fig. 
2c (g — 0.05eV), the shape of the zero-energy Landau 
band is similar to the clean case and T broadened by 
the disorder is as small as 3.6meV. In this case, the 
dominant factor on the subband shape is still the artifi- 
cial cut-off energy. However, the broadening width T of 
zero-energy Landau subband becomes large for g > O.leV 
when the disorder dominate the band broadening and the 
artificial cut-off energy is irrelevant. Numerical results 
show that the broadening shape of zero-energy Landau 
subband changes also from the Lorentizian shape to the 
Gaussian one. Figs. 2d-2e show that DOS can be fitted 
well by the Gaussian function, suggesting DOS of zero- 
energy Landau subband is Gaussian in the presence of a 
white noise potential^, [2f|. This finding is very simi- 
lar to the Gaussian-like shape for the lowest LL for 2D 
Schrodinger electrons @, ©, M, [H, E 0, 0, EE | . 
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FIG. 3: (Color online) Disorder-dependence of DOS for fixed 
B = 5T and B = 20T. (a) T vs. g for n = 0- and n = 
1— bands. It shows weak (or no) depednence of T on band 
index n. The linear relationship suggests V oc g. (b) DOS of 
zero-energy band for various g and at B = 20T. (c) 7) (0) vs. 
l/<? for the zero-energy band at B = 20T. 



FIG. 4: (Color online) Field-dependence of DOS for various 
fixed g. (a) r vs. \f~B for n = 0— and n = 1— bands. (? 
is fixed at 0.2eV and 0.3eV. The data shows weak (or no) 
depednence of F on band index n. The linear relationship 
suggests r oc y/B. (b) DOS of zero-energy band for various B 
and at a fixed g = 0.3eV. (c) D(0) vs. \[B for the zero-energy 
band for g — 0.3eV. 



As shown in Fig. 2b, Landau subbands are distorted 
and shifted toward zero energy direction in the pres- 
ence of white- noise disorder [2 7j except the zero-energy 
subband because disorder does not destroy the electron- 
hole symmetry. This disorder effect can be seen clearly 
from the noticeable move of the DOS peak positions of 
bands. This shift of DOS peaks can be qualitatively ex- 
plained within the Self-consistent non-crossing approxi- 
mation (SCNCA) According to SCNCA, DOS peak 
positions should move toward zero energy in the pres- 
ence of white-noise disorder due to the finite real part of 
the self energies, and the shift increases with g. How- 
ever, SCNCA fails to account for the shape of broadened 
bands. 

Fig. 3a shows the g— dependence of subband width T 
for the zero-energy subband (n = 0) for B = 5T and 20T. 
As expected, T of zero-energy Landau subband increases 
with g. The same dependence is obtained for the first 
(n = 1) Landau subband (squares). The data shows ap- 
proximately a linear dependence in g. The slop depends 
on the magnetic field. However, the width does not de- 
pends on the band index n as suggested by the complete 
overlap of T for n — (circles) and for n = 1 (squares) 
at B = 5T. Fig. 3b shows the DOS of zero-energy sub- 
band at B — 20T for various disorder g. The peak D(0) 
decreases with <?, and is proportional to the inverse of g 
as shown by the good linear fit in Fig. 3c. This is be- 
cause the total number of states in each Landau subband 
is a constant which is the total number of flux quanta in 
the system. Since T is proportional to g, to keep the 
area defined by the DOS curve and E-axis to a constant, 
the peak must be roughly follow 1/g dependence. More 
interestingly, the shape of zero-energy subband remains 
symmetric due to electron-hole symmetry, and may fol- 



low a true Gaussian distribution. 

Fig. 4 is the magnetic field dependence of DOS at a 
fixed g. Fig. 4a shows square-root field dependence of 
r for both n = 0— band and n = 1— band. The good 
overlap of T for the two subbands at g = 0.2eV suggests 
that, similar to the T — g dependence, T — B curves are 
universal (independent on band indices). Just like DOS-g 
dependence, Fig. 4b shows how DOS of zero-energy band 
change with magnetic field. It is clear that D(0) is pro- 
portional to square-root of B as shown in Fig. 4c. This 
is because the total number of states in one Landau sub- 
band is proportional to B, thus the product of the band 
peak (-D(Q)) and bandwidth T should be proportional to 

Vb. 

In order to understand our numerical findings quan- 
titatively, we generalize the approach of Wegner[8( for 
Schrodinger electrons to the Dirac electrons to obtain an 
analytic expression of DOS for the zero-energy band since 
it is only valid for small E. In the high magnetic field 
regime, the Green's function of unperturbed system for 
the zero-energy band is 



G (r,r',E) =< r\- 



1 



H 



\r' >= C(r,r') 



1 



E-En 



■(3) 



with 



C(r, r') = poexpi-(Ct - 2£'*£ + tO/^hf) 




1 



where £ = x + iy and £' = x' + iy' are the complex vari- 
ables, po = (2ttII)^ 1 is the degeneracy per area in each 
LL for the clean graphene. If the intervalley scattering 
between K and K' is ignored, the Green's function of 
the zero-energy LL has only one non-zero term, which 
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simplifies the problem a lot. In this case, the shape of 
zero-energy LL can be directly mapped to the problem of 
conventional 2DEG for the lowest LLQ. Therefore, the 
shape of the zero-energy band has a Gaussian-like form: 



/cn Po d . -I/ 2 



ry/E 2 /Po9 2 



7T 



'0 



2 dr,). (4) 



This Gaussian-like form is little bit flater than the true 
Gaussian distribution. Although our numerical date can 
also fit this Gaussian-like form, but they fit the Gaussian 
function better. This may be due to the fact that in- 
tervalley scattering is neglected in our derivation of the 
Gaussian-like form while our numerical calculation in- 
clude all effects from both disorders and magnetic field. 

The shape of Landau bands is important to many phys- 
ical quantities. As an example, we would like to evalu- 
ate the magneto-conductivity of a band whose DOS is 
a Gaussian function. In the Drude-Boltzemann approxi- 
mation, the Kubo formula yields the diagonal magneto- 
conductivity for the Fermi energy at the Dirac point [30j|: 



ie 2 te- 



ll 1 



(5) 



where t = » ^ and lu c — Vp/l c is cyclotron frequency. 
Although Eq. (4) works well only in the high magnetic 
field (t > 1). Noted that t is independent on B as 
r oc gVB, thus magneto-conductivity does not depend 
on the magnetic field, and is only disorder dependent. 
This behavior reflects that a xx should tend to saturate 
in high magnetic field, which is consistent with the ex- 
perimental observations of Ref. 31 1 . 

In conclusion, the shape of disorder broadened zero- 
energy Landau band can be described well by a Gaussian 
function. For white-noise random potential, band width 
is proportional to both the square-root of magnetic field 
and variance of random potential. The perturbation the- 
ory is used to justify the conclusions. The importance 
of the shape of Landau bands is demonstrated by the 
magneto-conductivity. The results can shield some lights 
on the properties of graphene in magnetic field. 
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